Vortex-line percolation in the three-dimensional 
complex IV^I^ model 
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In discussing the phase transition of the three-dimensional complex \ip\'^ theory, we study the ge- 
ometrically defined vortex-loop network as well as the magnetic properties of the system in the 
vicinity of the critical point. Using high-precision Monte Carlo techniques we investigate if both 
of them exhibit the same critical behavior leading to the same critical exponents and hence to a 
consistent description of the phase transition. Different percolation observables are taken into ac- 
count and compared with each other. We find that different connectivity definitions for constructing 
the vortex- loop network lead to different results in the thermodynamic limit, and the percolation 
thresholds do not coincide with the thermodynamic phase transition point. 
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I. INTRODUCTION 

Substantial progress in the understanding of the na- 
ture of phase transitions driven by topological excitations 
has been achieved in the beginning of the 1970's, when 
Berezinskiii and Kosterlitz and Thouless^- published their 
seminal papers on the two-dimensional XY model, in- 
volving the unbinding of point like vortices when the 
temperature exceeds a critical value. A few years later in 
1977, Banks, Myerson and Koguli^ showed that the Vil- 
lain model in three dimensions, a particular spin model 
with global 0{2) symmetry due to the 27r-periodicity 
in the Hamiltonian, can be represented by an equiva- 
lent defect model with long-range Biot-Savart-like inter- 
actions, where the spin configurations are integer valued 
and sourceless. These configurations can be interpreted 
as line like excitations forming closed networks which can 
be identified with the vortex loops of the original theory. 
At the transition point, where the broken 0(2) symme- 
try in the low-temperature phase is restored, loops of in- 
finite length become important which provides the basis 
for attempting a percolational treatment,'* So the ques- 
tion arises whether the percolational threshold coincides 
with the thermodynamic critical point, or under which 
conditions such a coincidence can be established.^ 

Percolational studies of spin clusters in the Ising model 
showed that one has to handle this approach carefully. 
It only works, if one uses a proper stochastic defini- 
tion of clusters iSiSiS Such so-called Fortuin-Kasteleyn 
clusters of spins can be obtained from the geometrical 
spin clusters, which consist of nearest-neighbor sites with 
their spin variables in the same state, by laying bonds 
with a certain probability between the nearest neighbors. 
The resulting Fortuin-Kasteleyn clusters are in general 
smaller than the geometrical ones and also more loosely 
connected. It is well known that depending on the clus- 
ter type considered one may find different sets of critical 
exponents and even different percolation thresholds, so a 
careful treatment is required. 

In three-dimensional, globally 0(2) symmetric theories 



the percolating objects are vortex lines forming closed 
networks. The question we want to address in this pa- 
per is: Is there a similar clue in the case of vortex net- 
works as for spin clusters, or do they display different 
features? Therefore we connect the obtained vortex line 
elements to closed loops, which are geometrically defined 
objects. When a branching point, where n > 2 junc- 
tions are encountered, is reached, a decision on how to 
continue has to be made. This step involves a certain 
ambiguity. We want to investigate the influence of the 
probability of treating such a branching point as a knot. 
This work concentrates on the three-dimensional complex 
Ginzburg-Landau model, the field theoretical representa- 
tive of the 0(2) universality class. 

The rest of the paper is organized as follows. In SecHll 
we give the definition of the model and introduce the 
observables. The results of our Monte Carlo simulations 
are presented in Sec. lIIII and concluding remarks and an 
outlook to future work can be found in Sec. IIVI 



II. MODEL AND OBSERVABLES 

The standard complex or two-component Ginzburg- 
Landau theory is defined by the Hamiltonian 



7>0, (1) 



where = "(pxir) + i^'yi'^) = |'0(7')|e*"^^''^ is a com- 

plex field, and a, b and 7 are temperature independent 
coefficients derived from a microscopic model. In order 
to carry out Monte Carlo simulations we put the model 
on a c?-dimensional hypercubic lattice with spacing 
a. Adopting the n otation of Ref. fiol we introduce scaled 
variables tp = tp/ ^ (|a|/6) and u = r/^, where ^ = 7/|a| 
is the mean-field correlation length at zero temperature. 
This leads to the normalized lattice Hamiltonian 
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where fi denotes the unit vectors along the d coordinate 
axes, N = L'^ is the total number of sites, and an unim- 
portant constant term has been dropped. The param- 
eter Vq merely sets the temperature scale and can thus 
be absorbed in the definition of the reduced temperature 
f = T/Vo- 

After these rescalings and omitting the tilde on a, 
and T for notational simplicity in the rest of the paper, 
the partition function Z considered in the simulations is 
given by 



Z : 
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where (3 — 1/T denotes the inverse temperature and 
J Dtp Dip = J D Retp D Im-0 stands short for integrating 
over all possible complex field configurations. 

In the limit of a large parameter tr, it is easy to read 
off from Eq. Q that the modulus of the field is squeezed 
onto unity such that the XY model limit is approached 
with its well-known continuous phase transition in three 
dimensions at (3c ~ 0.45?^ 

In order to characterize the transition we have mea- 
sured in our simulations to be described in detail in the 
next section among other quantities the energy (H), the 
specific heat = ((77^) — (H)'^) /N, and the mean-square 

amplitude = E^=i(IV'«P). In order to de- 

termine the critical temperature, the helicity modulus, 

1 ^ 
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and the Binder cumulant U = (| A/|^)/(|A/|^}^ were also 

computed, where M — + iMy = X]^=i V'n is the 
magnetization of a given configuration. 

The main focus in this paper is on the properties of the 
geometrically defined vortex-loop network. The standard 
procedure to calculate the vorticity on each plaquette is 
by considering the quantity 
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where (pi, . . . ,(f)4 are the phases at the corners of a pla- 
quette labeled, say, according to the right-hand rule, 
and [a\2TT stands for a modulo 27r: [a]27i- = a + 27m, 
with n an integer such that a + 2Tm G (— 7r,7r], hence 
m — ni2 -l- 7123 + i^3i + ^41- If TO 7^ 0, there exists a topo- 
logical charge which is assigned to the object dual to 
the given plaquette, i.e., the (oriented) line elements 
which combine to form closed networks ( "vortex loops" ) . 



FIG. 1: If two (or three) vortex lines pass through one 
cell, the vortex tracing algorithm must decide how to connect 
them, and this leads to an ambiguity in the length distribu- 
tion. Left: Connecting all line elements (forming a knot). 
Right: Connections are made stochastically. 



With this definition, the vortex "currents" can take 
three values: 0, ±1 (the values ±2 have a negligible prob- 
ability and higher values are impossible). The quantity 
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serves as a measure of the vortex-line density. 

In order to study percolation observables we connect 
the obtained vortex line elements to closed loops, which 
are geometrically defined objects. Following a single line, 
there is evidently no difficulty, but when a branching 
point, where n > 2 jmictions arc encountered, is reached, 
a decision on how to continue has to be made. This step 
involves a certain ambiguity. If we connect all in- and 
out-going line elements, knots will be formed. Another 
choice is to join only one incoming with one outgoing line 
element, with the outgoing direction chosen randomly. 
These two possibilities are shown in Fig. ^ We will em- 
ploy two "connectivity" definitions here: 

• "Maximal" rule: At all branching points, we con- 
nect all line elements, such that the maximal loop 
length is achieved. That means each branching 
point is treated as a knot. 

• "Stochastic" rule: At a branching point where 
n > 2 junctions are encountered, we draw a uni- 
formly distributed random number G (0, 1] and if 
this number is smaller than the connectivity param- 
eter c we identify this branching point as a knot of 
the loop, i.e., only with probability < c < 1 a 
branching point is treated as a knot. In this way 
we can systematically interpolate between the max- 
imal rule for c = 1 and the case c = 0, which corre- 
sponds to the procedure most commonly followed 
in the literature^ 

We can thus extract from each lattice configuration a 
set of vortex loops, which have been glued together by 
one of the connectivity definitions above. In Fig. |21 we 
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FIG. 2: (Color online) Vortex-loop networks for the "stochastic" definition with c = 0.4 (left) and the "maximal" definition 
with c — 1 (right). Both networks are generated from the same L = 8 lattice field configuration at the (thermodynamically) 
critical coupling /3c = 0.780 08. The different loops are distinguished by the color coding. 



show two possible vortex-loop networks for c = 0.4 and 
c = 1 generated out of the same lattice configuration. 

For each loop in the network, we measure the following 
observables: 

• "Mass", Oinass: The "mass" of a vortex loop is the 
number of line elements *l^^n of the loop, i.e., sim- 
ply its length Zioop normalized by the volume 

O 

mass — 

hoop/N . (8) 

By summing over all loops of a configuration we 
recover of course the vortex density Q, 

Oniass = V. (9) 

loops 

For the percolation analysis the mass of the longest 
loop Of^l^g in each vortex network is recorded, 
which usually serves as a measure of the percolation 
strength (behaving similarly to a magnetization)^ 

• "Volume", Ovoi: For each vortex loop, first the 
smallest rectangular box is determined that con- 
tains the whole loop. This value is then normalized 
by the volume of the lattice. A vortex loop spread 
over an extent l^, ly, and Iz thus results in 

OVOI - [Ix XlyX U)/N . (10) 

For each lattice configuration, we record the max- 
imal "volume" O^^, which may be taken as an 
alternative definition of the percolation strength. 

• "Extent" of a vortex loop in 1, 2, or 3 dimensions, 
OiY),02i), and Osr,: This means simply to project 
the loop onto the three axes and record whether 
the projection covers the whole axis, or to be more 



concrete, whether one finds a vortex-line element 
of the loop in all planes perpendicular to the eyed 
axis. If there is a loop fulfilling this requirement, 
then this loop is percolating and we record 1 in the 
time series of measurements; if not, a value of is 
stored. This quantity can thus be interpreted as 
percolation probability^! which (behaving similarly 
to a Binder parameter) is a convenient quantity for 
locating the percolation threshold f3p. 

• "Susceptibilities," Xi- For the vortex- line density v 
and any of the observables Oi defined above (i = 
"mass," "vol," "ID," . . . ), one can use its variance 
to define the associated susceptibility, 

X^ = Ni{0^)^(0,r) , (11) 

which is expected to signal critical fluctuations. 

• "Line tension," 6: On general grounds the loop- 
length distribution P{hoop) is expected to have the 
following form:^^ 

^aioop)^Cpexp(-?ioop0) , (12) 

where the Fisher exponent r is given in terms of 
the fractal dimension D of the loops by 

r=^ + l. (13) 

For a three-dimensional (d = 3) (nonintcracting) 
Brownian random walk with D = 2 this leads 
to T = 5/2, while r > 5/2 for self-avoiding and 
r < 5/2 for self-seeking lines, respectively. The 
parameter 9 is the line tension which vanishes ac- 
cording to^ 

^ = l/3-/3p^^ (14) 
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where /3p is the percolation threshold of the ran- 
dom walk and 79 = l/crg the second independent 
percolation exponent 1^ 

III. SIMULATION AND RESULTS 

Let us now turn to the description of the Monte Carlo 
update procedures used by us. We employed the single- 
cluster algorithm^^ to update the direction of the fieldfi^ 
similar to simulations of the XY spin model. The mod- 
ulus of V' is updated with a Metropolis algorithmii^ii& 
Here some care is necessa ry t o treat the measure in 
Eq. @ properly (see Ref. llTll . One sweep consisted 
of N spin flips with the Metropolis algorithm and Ngc 
single-cluster updates. For all simulations the num- 
ber of cluster updates was chosen roughly proportional 
to the linear lattice size, A'sc ~ L, a standard choice 
for three-dimensional systems as suggested by a simple 
finite-size scaling (FSS) argument. We performed sim- 
ulations for lattices with linear lattice size L = 6 - 
10,12,14,16,18,20,22,24,26,28,32,36, and 40, respec- 
tively, subject to periodic boundary conditions. After 
an initial equilibration time of 20 000 sweeps we took 
about 100 000 measurements, with ten sweeps between 
the measurements. All error bars are computed with the 
Jackknife method.^^ 

In order to be able to compare standard, thermody- 
namically obtained results (working directly with the 
original field variables) with the percolative treatment 
of the geometrically defined vortex-loop networks con- 
sidered here, we used the same value for the parameter 
a = 1.5 as in Ref.figl for which we determined by means 
of standard FSS analyses a critical coupling of 

/3c = 0.780 08(4) . (15) 

Focussing here on the vortex loops, we performed new 
simulations at this thermodynamically determined criti- 
cal value, /3 = 0.780 08, as well as additional simulations 
at /3 = 0.79, 0.80, and 0.81. The latter /3 values were 
necessary because of the spreading of the pseudocritical 
points of the vortex loop related quantities. As previ- 
ously we recorded the time series of the energy i/, the 
magnetization M, the mean modulus |-0|, and the mean- 
square amplitude I V'P, as well as the helicity modulus 
and the vortex-line density v. In the present simulations, 
however, we saved in addition also the field configura- 
tions in each measurement. This enabled us to perform 
the time consuming analyses of the vortex-loop networks 
after finishing the simulations and thus to systematically 
vary the connectivity parameter c of the knots. 

The FSS ansatz for the pseudocritical inverse temper- 
atures l3i{L), defined as the points where the various Xi 
obtain their maxima, is taken as usual as 

A {L) = /3.,c + ciL-i/"^ + + . . . , (16) 

where (3i^c denotes the infinite- volume limit, and 1/ and uj 
are the correlation length and confluent correction crit- 



ical exponents, respectively. Here we have deliberately 
retained the subscript i on (3i^c- 

Let us start with the susceptibility Xv of the vortex- 
line density. Note that this quantity, while also being ex- 
pressed entirely in terms of vortex elements, plays a spe- 
cial role in that it is locally defined, i.e., does not require 
a decomposition into individual vortex loops (which, in 
fact, is the time-consuming part of the vortex- network 
analysis). Assuming the XY model values for v and u 
compiled in Table HI which are taken from Refs. and 
I20I and fitting only the coefficients (3i^c and Cj, we arrive 
at the estimate 

/3,,c = 0.7797(14) (17) 

with a goodness-of-fit parameter Q = 0.20. This value is 
perfectly consistent with the previously obtained "ther- 
modynamic" result H15() . derived from FSS of the mag- 
netic susceptibility and various (logarithmic) derivatives 
of the magnetization. On the basis of this result, it would 
be indeed tempting to conclude that the phase transition 
in the three-dimensional complex Ginzburg-Landau field 
theory can be explained in terms of vortex-line prolifer- 
ationmii^ As pointed out above, however, the vortex- 
line density v does not depend on the connectivity of the 
vortex network and therefore does not probe its percola- 
tion properties. In fact, v behaves similar to the energy 
and the associated susceptibility Xv similar to the spe- 
cific heat, so that the good agreement between Eqs. H15(l 
and H17|l is a priori to be expected. 

To develop a purely geometric picture of the mecha- 
nism governing this transition, one should thus be more 
ambitious and also consider the various quantities Oi in- 
troduced above that focus on the percolative properties of 
the vortex-loop network. As an example for the various 
susceptibilities considered, we show in Fig. O the suscep- 
tibility X3D of OsD for c = and c = 1. The resulting 
scaling behavior of the maxima locations /33d(^) is de- 
picted in Fig. 21 where the lines indicate fits according 
to Eq. (|16|1 with exponents fixed again according to Ta- 
ble d We obtain /Jgo.c = 0.7824(1) with xVdof = 1.14 
(Q = 0.32, L>8) for c = and /Jgo.c = 0.8042(4) with 
xVdof = 0.75 (Q = 0.58, L > 20) for c = 1. While for 
the "stochastic" rule with c — the infinite- volume limit 
of (33b{L) is at least close to (3c, it is clearly significantly 
larger than Pc for the fully knotted vortex networks with 
c 1. 

By repeating the fits for all vortex-network observables 
and the parameter c between and 1 in steps of 0.1, we 
find the results collected in Tables ITU and UTTl To check 
the stability of the fit results we performed fits with dif- 
ferent lower bounds of the fit range Lmin, while the upper 
bound was always our largest lattice size L = 40. For all 
observables, except for Osd, we found a weak dependence 
ofPi^c on the fit range. For all five observables we see that 
the location of the infinite- volume limit Pi^c does depend 
on the connectivity parameter c used in constructing the 
vortex loops in a statistically significant way. With de- 
creasing c, the infinite- volume extrapolations come closer 
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TABLE I: The critical exponents of the 3D XY model universality class as reported in Ref. and the correction-to-scaling 
exponent uj of Ref. IT^ 



a (3 'y S rj u Lu 

-0.0146(8) 0.3485(2) 1.3177(5) 4.780(2) 0.0380(4) 0.67155(27) 0.79(2) 



toward the thermodynamical critical value H15|l , but even 
for c = theyclearly do not coincide. 

As in Ref. we found that the percolation points /3i.c 
of Oi satisfy some inequalities. Because each lattice cube 
has three plaquettes, Ovoi > C'mass/S, and it is plausible 
that (Oid) > {02d) > (Osd). The first relation implies 

mass.c • 

(18) 

Our results collected in Table ITll are consistent with this 
inequality. In addition to this inequality the authors of 
Ref.i also conjectured that Pvo\,c = /33D,c = /3id,c- Our 
numerical data show that /3voi,c ~ Piu.c, but the other 
percolation points satisfy only the following inequalities: 

/3vol,c ~ /?lD,c > /32D,c > /33D,c , (19) 

cf. Table lllil The reason for this are possibly different 
corrections to scaling for the different observables. In 
the infinite- volume limit all definitions should lead to the 
same critical point. 

These findings are reminiscent of the percolation be- 
havior of, say, Ising (minority) spin droplets of like spins 
which are known to percolate in three dimensions already 
below the transition temperature, i.e., (3p > Pc as for 
the vortex-loop observables. Only by breaking bonds be- 
tween like spins with a certain temperature dependent 
probability (= exp(— 2/3)), one can tune the thus de- 
fined Fortuin-Kasteleyn (FK) clusters to percolate at (3c- 
With any other non-FK probability < pb < p^^ for 
breaking bonds between like spins it is conceivable that 
the associated percolation point would be located some- 
where between /3p of the geometrical droplets and the 
thermodynamical (or, equivalently, FK) critical point (3c 
(for pb > the percolation transition may even vanish 
altogether). By analogy, our connectivity parameter c 
seems to play a similar role for the vortex-loop network 
as pb for the spin droplets. However, due to the miss- 
ing analog to the FK representation of the Ising model, 
in the present case of the vortex-loop network, it is not 
easy to guess a suitable temperature dependence of the 
parameter c and we hence eluded to using a systematic 
variation of c in small constant increments. The other 
important difference to the case of Ising droplets is of 
course the long-range interaction between vortex-line el- 
ements which certainly puts the sketched analogy to Ising 
droplets on quite an uncertain and speculative footing. 

With these remarks in mind we nevertheless performed 
tests whether at least for c = the critical behavior of 
the vortex-loop network may consistently be described 
by the three-dimensional XY model universality class. 



As an example for a quantity that is a priori expected to 
behave as a percolation probability we picked again the 
quantity O3D for which the susceptibility was already 
shown in Fig. |3| As is demonstrated in Fig. [Sfa) for the 
case c = 0, by plotting the raw data of O3B as a func- 
tion of (3 for the various lattice sizes, one obtains a clear 
crossing point so that the interpretation of O3D as perco- 
lation probability is nicely confirmed. To test the scaling 
behavior we rescaled the raw data in the FSS master plot 
shown in Fig.|SJb), where the critical exponent v has the 
XY model value given in Table |l| and /3c(C3d) = 0.7842 
was independently determined by optimizing the data 
collapse, i.e., virtually this is the location of the cross- 
ing point in Fig.|SJa). The collapse turns out to be quite 
sharp which we explicitly judged by comparison with sim- 
ilar plots for standard bond and site percolation (using 
there the proper percolation exponent, of course). For 
c > we found also a sharp data collapse, but for a 
monotonically increasing exponent v, which is for large c 
values compatible with the percolation critical exponent 
v — 0.8765(16) on a three-dimensional simple cubic lat- 
ticeJS One should keep in mind, however, that neither 
/SsD.c as extrapolated from the susceptibility peaks nor 
the estimate obtained from the crossing point in Fig.[SJa) 
is compatible with Pc- 

Next we looked at Omass which a priori is expected to 
behave like a percolation strength, that is similarly to the 
magnetization with an inverted /? axis. The plot of the 
raw data for c = as a function of (3 in Fig. IHIa) indeed 
seems to confirm this expectation. To test the scaling 
properties we show in Fig. IH^b) the corresponding FSS 
master plot, where the critical exponents v and (3 are 
again fixed to their XY model values (cf. Table UJ and 
/3c(Cmass) = 0.782 75 was determined by optimizing the 
data collapse. Also this collapse is comparatively sharp. 
Even though the thus obtained value for /9c(Cmass) is con- 
sistent within error bars with the FSS value in Table UTI 
obtained from the susceptibility maxima locations (but 
even further away from (3c), we found a visible spread of 
the rescaled curves when the latter value was used and 
kept fixed. Similarly, assuming both XY model expo- 
nents and /3c (Omass) = /3c does not produce a satisfac- 
tory data collapse. Thus for both observables, O3D and 
Cmass, we obtain nice FSS scaling plots at c = compat- 
ible with XY model critical exponents, but a "wrong" 
critical coupling. 

Surprisingly, when using c ^ for constructing the 
vortex loops, Omass shows a completely different behav- 
ior. As example, we show in Fig. our data for the case 
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c = 0.4. Already by looking at the raw data in Fig. |7{a), 
it is obvious that, for c 7^ 0, the mass of vortex loops no 
longer behaves as a percolation strength (i.e., magneti- 
zation); rather it resembles pretty much the percolation 
probability O3D • From the crossing point of the curves we 
get /3(Omass) — 0.786 46. Using this value we get in the 
FSS master plot shown in Fig. [7][h) a nice data collapse 
for 1/ = 0.98, but in contrast to the c = case now only 
when the y axis is not rescaled. We repeated this analysis 
for all values c > and found a monotonically increas- 
ing exponent ^ from v « 0.74 for c = 0.1 to « 1.82 
for c = 1.0, which appears quite nonsensical. A precise 
determination of the critical exponents as a function of c 
was not our aim here and anyway, due to the rather small 
lattice sizes studied, also not feasible. Still, this strange 
behavior clearly calls for an explanation. 

The variance definition l|ll() of the susceptibilities stud- 
ied so far is quite unusual in percolation theory. We have 
therefore also investigated the standard percolation def- 
inition for the average loop size x/ioop (^^ ^^en at a given 
link of the lattice) which is expected to scale as the vari- 
ances defined above. In terms of the loop-length distri- 
bution P(Zioop) it is given as4 

^ _ Sloop ^oop^(^loop) ^^^^ 
Sloop ^loop-P(^loop) 

where the prime on the sum is to indicate that we dis- 
card in each measurement the percolating loop according 
to the criterion Osd- For this observable we also found 
a clear displacement between the maxima for different 
values of the connectivity parameter c and the thermo- 
dynamic transition point. Unfortunately, the reweighting 
range for c = was too narrow for this observable to al- 
low more detailed analyses, see Fig. |S1 For c > the 
location of the pseudocritical points of the average loop 
size behave similar to the susceptibilities as defined in 
Eq. (|ll|l and lead to slightly higher f3c values than the 
thermodynamic one. 

Finally, in Fig. IHIa) we show the loop-length distribu- 
tion P(^ioop) (without the largest loop) as a function of 
the loop length /loop for the "stochastic" rule (c = 0) 
for various temperatures and L = 40. From Eq. H12|) 
one expects that the decay changes from exponential to 
algebraic at (3c, because of the vanishing of the vortex- 
line tension 9. Also from this analysis we found that the 
percolation transition takes place at a slightly higher /? 
value than the thermodynamic one. We performed fits 
according to P(/ioop) ~ ^loop ^^'^ found the best results 
for P = 0.79 where r « 2.2 < 5/2, see TablelTVI We want 



to note that this /? value is for our largest lattice and we 
only examined the loop-length distributions at the /? val- 
ues used for the simulation. To determine the percolation 
transition and also the line tension with the help of the 
loop-length distributions one would need a finer temper- 
ature spacing. At Pc of the thermodynamic transition we 
looked at the change of the decay of the distributions as 
a function of the connectivity parameter c for L = AO, 
see Fig. EJ^b). Here we found the best agreement with 
an algebraic decay for c = 0.1 with t = 2.348(1) and 
X^/dof = 2.9. From this observation and the fact that 
we found a pronounced peak for the largest loops in the 
distribution at c = 0.1 and no peak at c = 0.0, we con- 
clude that the line tension vanishes close to c = 0.0. 



IV. CONCLUSION AND OUTLOOK 

In this paper we have found for the three-dimensional 
complex Ginzburg-Landau field theory that the geomet- 
rically defined percolation transition of the vortex-loop 
network is close to the thermodynamic phase transition 
point, but does not coincide with it for any connectivity 
definition we have studied. Our results for the connec- 
tivity parameter c S [0, 1] extend the claim of Ref. for 
the three-dimensional XY spin model that neither the 
"maximal" (c — 1) nor the "stochastic" rule (c — 0) 
used for constructing macroscopic vortex loops does re- 
fiect the properties of the true phase transition in a strict 
sense. Nevertheless it may be possible to bring the per- 
colation transition closer to the thermodynamic one by 
using different vortex-loop network definitions, e.g., using 
a temperature-dependent or a size-dependent connectiv- 
ity parameter in analogy to the Fortuin-Kasteleyn defini- 
tion for spin clusters. To verify this presumption would 
be an interesting future project, but thereby one should 
first investigate the XY model which is much less CPU 
time-consuming. 
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TABLE II; FSS fits according to Eq. 1161 in tiie range Lmin to Lmax = 40 for the mass and volume order parameter of the 
vortex loops for various values of c (c = 1: maximally knotted), with exponents v and lo fixed according to Table Q] The 
thermodynamic transition is at /3c = 0.780 08(4) (Ref. 19). 













Cvol 




c 




-t^min 


xVdof 


Q 


/3c 




xVdof 


Q 


1.0 


0.8017(8) 


20 


1.53 


0.17 


0.8072(4) 


12 


1.86 


0.05 


0.9 


0.8016(6) 


16 


1.66 


0.11 


0.8048(3) 


14 


1.08 


0.37 


0.8 


0.7996(9) 


16 


1.55 


0.14 


0.8028(3) 


14 


1.40 


0.19 


0.7 


0.7976(5) 


14 


1.34 


0.21 


0.8005(4) 


16 


1.14 


0.33 


0.6 


0.7970(5) 


16 


1.10 


0.36 


0.7995(5) 


14 


0.57 


0.80 


0.5 


0.7954(6) 


18 


1.32 


0.24 


0.7963(3) 


16 


1.23 


0.28 


0.4 


0.7912(5) 


14 


0.94 


0.48 


0.7938(4) 


18 


0.82 


0.55 


0.3 


0.7887(9) 


16 


1.32 


0.23 


0.7924(4) 


18 


1.79 


0.09 


0.2 


0.7856(3) 


6 


1.62 


0.07 


0.7907(4) 


12 


0.28 


0.98 


0.1 


0.7834(4) 


10 


0.37 


0.96 


0.7875(2) 


12 


0.94 


0.49 


0.0 


0.7811(39) 


10 


1.26 


0.24 


0.7834(3) 


16 


0.92 


0.48 



TABLE III: Same as in Table UTI for the vortex parameters Oid, C2D, and Osd- 










ID 









2D 






C»3D 






c 


/3c 


^min 


xVdof 


Q 


/3c 




xVdof 


Q 


/3c 


^min 


xVdof 


Q 


1.0 


0.8066(5) 


16 


0.68 


0.68 


0.8051(3) 


20 


0.86 


0.50 


0.8042(4) 


20 


0.75 


0.58 


0.9 


0.8048(3) 


16 


0.68 


0.69 


0.8032(3) 


20 


0.72 


0.61 


0.8027(3) 


18 


0.37 


0.89 


0.8 


0.8030(2) 


16 


1.04 


0.40 


0.8027(4) 


14 


0.74 


0.66 


0.8009(2) 


18 


0.54 


0.77 


0.7 


0.8011(2) 


16 


0.90 


0.50 


0.7999(3) 


18 


1.29 


0.26 


0.7988(2) 


18 


0.61 


0.72 


0.6 


0.7992(3) 


16 


1.67 


0.11 


0.7976(2) 


18 


0.44 


0.85 


0.7968(2) 


18 


1.20 


0.30 


0.5 


0.7966(3) 


16 


0.94 


0.47 


0.7953(2) 


18 


0.42 


0.86 


0.7953(2) 


18 


0.72 


0.67 


0.4 


0.7938(3) 


20 


1.18 


0.32 


0.7940(2) 


14 


1.10 


0.36 


0.7928(2) 


12 


0.35 


0.96 


0.3 


0.7919(2) 


18 


1.11 


0.35 


0.7915(3) 


16 


0.66 


0.71 


0.79037(5) 


6 


0.54 


0.91 


0.2 


0.7907(2) 


12 


1.10 


0.37 


0.7890(2) 


16 


0.68 


0.69 


0.7881(1) 


6 


0.92 


0.53 


0.1 


0.7868(3) 


18 


1.38 


0.21 


0.7861(2) 


18 


0.49 


0.81 


0.7857(2) 


6 


0.60 


0.86 


0.0 


0.7837(2) 


12 


1.84 


0.08 


0.7824(5) 


18 


0.28 


0.94 


0.7824(1) 


8 


1.14 


0.32 



TABLE IV: Results for the Fisher exponent r for the "stochastic" definition (c — 0) at three inverse temperatures (3 = 1/T, 
assuming an algebraic decay of the loop-length distribution close to criticality. 



/3 


fit 


range 


r 


xVdof 


0.78008 


20 


-500 


2.261(1) 


3.1 


0.79 


20 


-800 


2.201(1) 


1.4 


0.80 


20 


-500 


2.275(1) 


26 
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5 10 15 20 25 30 35 40 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

FIG. 4: Location of the percolation thresholds determined from the maximum of susceptibility of O30 for c — and c = 1 as 
a function of (a) L and (b) L~^^" , respectively. The lines indicate fits according to Eq. II6II with v and lj fixed according to 
TableU] The horizontal dashed line shows the thermodynamically determined critical coupling /3c ~ 0.780 08(4). 
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FIG. 6: (Color online) (a) Mass of vortex loops as a function of inverse temperature /3 — 1/T for c = 0. (b) Rescaled data 
assuming 3D XY model critical exponents (cf. TableHJ and adjusting /3c(C'mass) = 0.782 75 for the best data collapse. 
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FIG. 8: (Color online) Average loop size Xiioop ^ ^ function of inverse temperature f3 = 1/T for c — 0, 0.2, and 1 (from left 
to right). 
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FIG. 9: (Color online) Loop-length distribution P{lioop) for our largest lattice (L = 40) as a function of the loop length, (a) 
The "stochastic" rule (c = 0) for various temperatures, (b) Behavior at the thermodynamic critical point (3 — 0.780 08 for 
various values of the connectivity parameter c. At c ~ 0.1 the decay changes from exponential to algebraic implying that the 
vortex- line tension 6 vanishes. 



